Approach to asymptotic behavior in the dynamics of the trapping reaction 
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We consider the trapping reaction A + B — > B in space dimension d = 1, where the A and B 
particles have diffusion constants Da, Db respectively. We calculate the probability, Q(t), that 
a given A particle has not yet reacted at time t. Exploiting a recent formulation in which the B 
particles are eliminated from the problem we find, for t — > oo, Q(t) ~ exp[— (4/y / 7r)(p 2 _Dst) 1 ' /2 — 
(Cp 2 DAt) 1 ^ 3 + •••], where p is the density of B particles and C oc Da/Db for Da/Db <C 1. 
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I. INTRODUCTION 

The trapping reaction, A + B — > £?, where A and B 
are diffusing particles with diffusion constants Da and 
Db respectively, has seen a resurgence of interest in the 
last few years 0, 0, 0, S IE @ • The problem is related to 
the two-species annihilation reaction 0, A + B — > 0, 
in the case where the initial particle densities satisfy 
Pa(0) < Pb(0). In the long-time limit pA(t) -C Ab(*), 
and to study the asymptotic decay of the A-particle den- 
sity it suffices to study the survival probability, Q(t), of 
a single A particle diffusing in a sea of B particles. In 
this limit, the two-species annihilation reaction is equiv- 
alent to the trapping reaction. As the initial condition, 
the -B-particles are taken to be randomly distributed in 
space with uniform density p, i.e. the probability to find 
a B-particle in a volume element dV is pdV. 

The problem of determining the asymptotic behavior 
of Q(t) has only recently been solved [2-Lj- The solution 
is limited to space dimensions d < 2. The asymptotic 
form of the solution has been known for some time Q . 
The A-particle survival probability is given by 




-x d t d/2 ), 

-\ 2 t/lnt), 
-A d t), 



d < 2 
d= 2 
d > 2 



(1) 



but the values of the constants were not known. In 
refs. 0, 0] , these were determined for d < 2 by deriving 
upper and lower bounds on Q(t) and showing that they 
pinch as t — > oo. The restriction to d < 2 is a consequence 
of the recurrence properties of random walks. For d < 
2, all particles can be treated as zero-size points. As a 
consequence, the combinations p(DBt) d ^ 2 and Da/Db 
are the only dimensionless quantities in the theory. A 
key result of 0, LI is that the coefficients Ad in Eq. (JTJ 
are independent of Da for d < 2, and therefore have the 
form Ad = a,d pD d g 2 , where the quantities ad are pure 
numbers. The result of ^ s 



a d 



(4-7r) d / 2 (2/nd) sm(nd/2), d<2 
4tt, d = 2 



(2) 



For d — 2, the A-particle has to be given a non-zero size, 
but this size only enters as a cut-off in the logarithm in 
Eq. (UJ. 



Although these results are asymptotically exact, nu- 
merical studies 0, 13 show that the asymptotic regime 
is not reached, either in d = 1 or d = 2, after any time 
accessible to simulation, even using sophisticated simu- 
lation methods that allow values of Q(t) as small as 
1CP 70 to be achieved. It is natural, therefore, to try to 
calculate the leading correction to the asymptotic behav- 
ior. This is the principal aim of this paper. It should be 
stated at the outset that this is a very challenging prob- 
lem, and we have obtained concrete results only in the 
limit D A < D B - 

The starting point of our approach is an equation, first 
presented in that implicitly determines the survival 
probability of the A-particle for a given A-particle tra- 
jectory, averaged over all the B-particle trajectories and 
-B-particle initial conditions. In this equation, therefore, 
the -B-particles have been eliminated from the problem. 
Their presence is felt only through their density, p, and 
their diffusion constant, Db, which enters the equation 
through the _B-particle diffusion propagator. The final 
step in the calculation is to average the survival proba- 
bility over the A-particle trajectories, weighted with the 
Wiener measure. 



In the present work we focus on computing the leading 
pre-asymptotic corrections to the known asymptotic be- 
havior. For simplicity we specialize to dimension d = 1, 
where the best numerical results are available. Including 
the new correction, the result we obtain is fully consistent 
with the available data. 



The paper is organized as follows. In section II we 
set up the necessary formalism, eliminating the B parti- 
cles from the problem. In order to make the paper self- 
contained, in section III we briefly review the calculation 
of the leading asymptotic behavior. The pre-asymptotic 
corrections are discussed in section IV, at leading nontriv- 
ial order in the ratio Da/Db- In the short Section V we 
discuss how one might go beyond this leading order, and 
propose a general form for the pre-asymptotic correction 
that interpolates between small and large Da/Db- The 
paper concludes with a brief discussion and summary of 
the results. 
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II. ELIMINATING THE B-PARTICLES 

The key step in deriving the 'fundamental equation' 
from which we start is a conceptual one. We treat 
the A and B particles as if they were non-interacting, 
and exploit the initial condition that each B particle is 
randomly located anywhere in the system. Consider, 
for simplicity, a system of finite length, L, containing 
N = pL B-particles (diffusion constant D B ), randomly 
distributed within it, and a single A particle (diffusion 
constant Da), initially located at the origin. Let z{t) be 
the A-particle's trajectory, and let P(x,t) be the proba- 
bility that a given B particle, starting at x, has met the A 
particle before time t. The average of this quantity over 
the initial position, x, is (1/L) J dxP{x,t) = R(t)/L, 
where R(t) is an implicit functional of z{t). The proba- 
bility that n distinct B particles have met the A par- 
ticle, averaged over their initial positions, is p n (t) — 
( N n )(R/L) n {l - R/L) N - n . Taking the limit N -> oo, 
L — > oo, with p = N/L and n held fixed, yields the Pois- 
son distribution, 



Pn = — r exp(-/i), 



(3) 



with p[z] = pR[z], the notation emphasizing that p is a 
functional of the A-particle trajectory z(t). The Poisson 
distribution implies that p[z] is the mean number of dis- 
tinct B particles that have met the ^4-particle up to time 
t. Its derivative, p(t), is the crossing rate: pit) dt is the 
probability that a B particle meets the A particle for the 
first time in the time interval (t, t + dt). 

The probability that the trajectory z(t) has survived, 
in the original interacting problem, is simply the proba- 
bility that there have been no crossings: 



Po[z] = exp(-p,[z]). 



(4) 



Finally, the survival probability Q(t) is obtained by aver- 
aging exp(— fi[z]) over all possible A-particle trajectories 
z(t) with the Wiener measure: 



Q(t) 



J Dz(t) exp 


-iBi/oMi) 2 -^] 


J Dz(t) exp 


1 f* J 1 dz\ 2 
AD A Jo U \dr) 





(5) 



where the path-integrals are restricted to functions z(t) 
satisfying z(0) = 0. The 'action functional' 



S[z] = 



1 



AD A 



+ p[z] 



(6) 



appearing in the numerator of Eq. defines an effective 
dynamics for surviving A-particle trajectories. 

The functional p[z] can be determined as follows. We 
calculate, in two ways, the probability density to find a B 
particle at the end-point of the trajectory, z(t), at time t. 
First, since the particles are treated as non-interacting, 
and the B particles start in a steady-state configuration 



of uniform density p, this probability density is simply p. 
It follows that 



p= / dt' p(t')G{z(t),t\z{t'),t') 



(7) 



where, on the right-hand side, p{t')dt' is the probability 
for a B particle to have its first encounter with the A in 
the time interval (t', t' + dt'), and 



G(z(t),t\z(t'),t') = 



1 



[4ttD b (t-t 1 )] V* 

[z(t)-z(t')f 



x exp 



AD B {t-t') 



(8) 



is the probability density for this £?-particle to subse- 
quently arrive at z{t) at time t. Eq. J7J is an implicit 
equation for the functional p[z] (noting that p(t = 0) = 
0, since no B particle can meet the A particle in zero 
time). The continuum version presented here was first 
derived in ref . . The remainder of the paper deals with 
some of the consequences of the 'fundamental equation' 
0. 



III. EXACT ASYMPTOTICS OF Q(t) 

Before dealing with the pre-asymptotic behavior, 
which is the main focus of this paper, we briefly review 
the derivation of the leading asymptotic behavior first 
presented in refs. 0- 0- [3] The argument is based on 
deriving upper and lower bounds for Q(t) and showing 
that these bounds pinch for t — > oo. 

First consider the case where the A particle is station- 
ary (i.e. Da = 0), the so-called 'target annihilation prob- 
lem', where z(t) = for all t. This problem can be easily 
solved by a variety of methods 0, . Here we treat it 
as a simple application of Eq. J7J). We put z = in Eq. 
(JHJ), and call the corresponding solution of J7J po(t). It 
satisfies the equation 



dt' p {t')[AvD B (t-t')]- l t 2 



(9) 



This equation is readily solved by taking the Laplace 
transform of both sides, since the right-hand side has 
the form of a convolution. Solving for po[s), and invert- 
ing the transform, the result is po(t) = Xit 1 / 2 , as in QJ, 

1 /2 __ 

with Ai = pDg a\ and a\ given by iJSJ. 

It is intuitively clear that Q(t) for the system with 
Da = (stationary A particle) provides an upper bound 
on Q(t) for a system with Da > 0. This bound can 
be proved rigorously Q using Eq. (Q. First we write 
p = po+pi in JJJ, and write, in a more compact notation, 



G(h,t 2 ) 



[iirDsitx-h)} 1 /- 



-[l-K{t x ,t 2 )], (10) 
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where we have taken t\ > t<i without loss of generality, 
and 

K{t x M) = l-exp{-[«(t 1 )-«(t 2 )] 2 /4I»fl(ti-t2)}. (11) 
Eq. J7J can then be rearranged in the form 

(12) 

where the full /i appears on the right-hand side. Taking 
Laplace transforms of both sides, solving for fli(s), and 
inverting the transform gives 



dti 



Wo (t-hyi* 



A(*2)if(tl,t2). (13) 



Eq. (|13fl is an implicit equation for fii{t) because the 
full fj, = fMo + fii appears on the right-hand side. Note, 
however, that K{t\,t2) > and (i > (because /j,(t) - 
the mean number of different B particles that have met 
the A particle up to time t — is clearly a non-decreasing 
function). Therefore n\[z] > for all paths z(t), with 
equality when z{t) = for all t. It follows that n[z\ > 
Ho(t) for all paths z{t), and 



Q{t) < Quit) = exp[-no(*)] • 



(14) 



An equivalent result on the lattice, in any space dimen- 
sion, was obtained in ref. |(J. 

This rigorous upper bound for Q{t) has the same 
asymptotic behavior as the rigorous lower bound derived 
m 01 > proving that the asymptotic form of Q{t) is 
the same as for the 'target'problem, where the A parti- 
cle is stationary. We sketch the derivation of the lower 
bound in d = 1. The results for other dimensions in the 
range d < 2 can be found in @, Q • To derive the d = 1 
result we can calculate the probability for a subset of his- 
tories for which the A particle survives. These are the 
histories in which (i) The region {—1/2, 1/2) initially con- 
tains no B particles; (ii) No B-particles enter this region 
up to time t; and (iii) the A particle stays within this 
region up to time t. The probabilities associated with 
(i)-(iii) are exp(— pi), exp[— fXo{t)], and (asymptotically) 
~ exp(— n 2 DAt/l 2 ) respectively. The product of these 
probabilities gives, for any 'box size' I, a lower bound on 
Qit). Optimizing this bound with respect to I, at any 
given time t, gives the best lower bound: 

Q{t) > Q L {t) ~ exp[-/i («) - 3(tt 2 p 2 -Da«/ 4 ) 1/3 ]- (15) 

Since /ito(i) ~ t 1 ^ 2 for d = 1, the second term in 115fl is 
subdominant for large t, and the two bounds agree. This 
proves that the target problem {Da = 0) gives the exact 
asymptotic form of Q{t) for any Da- Extensions of this 
result to general dimension d < 2 are given in Q . 

In the following section we discuss how one might cal- 
culate the leading pre- asymptotic correction to Q{t). 



IV. PRE-ASYMPTOTIC CORRECTIONS 

A. preliminaries 

The starting point of the calculation is Eq. JSJ, which 
is an exact expression for the survival probability, Q{t), 
of an A particle with diffusion constant Da moving in a 
sea of noninteracting B particles with diffusion constant 
Db- The functional fi[z] appearing in this equation is 
determined implicitly by the fundamental equation Q. 
It can be determined perturbatively as follows. Writing 
fj,[z] = fio[z] + Ml Hi where jjbQ = Aii 1 / 2 , we have shown 
that H\[z] is given by Eq. ffT^i where K{t\,t-z) is given 
by Eq. (|llfl . For small diffusion constant, Da, of the A 
particle, the typical excursion of the end point, z{t), of 
the A particle trajectory will be small, and the function 
K can be expanded in powers of z 2 : 

[z(ti)-*(f 2 )] 2 1 [z{h)-z{t^ 
[1,2> 4D B {h-t 2 ) 2[4D B {t 1 -t 2 )} 2+ "' 

This expansion is equivalent, as we shall see, to an expan- 
sion in powers of Da- Replacing K{t\,t2) in the right- 
hand side of Eq. Ijl3|l by the first term in the expansion 
lJT?)Jl. and /t(*2) by /iofe), where fi = Xit 1 / 2 , gives, to 
first order in Da, 



Ai 



dtx 



8ttD b J {t-H) 1 ' 2 
dU 



o 4 /2 (ii-< 2 ) 3 / 2 



Wi) ~ *{h)f 



(17) 



which is a quadratic functional of the ^4-particle trajec- 
tory z(t). Working at the quadratic level will be adequate 
to first order in Da- 



B. Survival probability 

Using the quadratic expression, Eq. (T7J|, for li\[z], we 
have the following expression for the survival probability: 



r>M ( \ +l/2^ I Dz{t)exp{~S[z}) 

Q{t) = exp -A^ / ) - (18) 

J Dz{t) exp(-S'o[2;J) 

where both path integrals are over all trajectories with 

initial points z{0) = 0, but with the end points, z{t), 

unconstrained. The 'action functionals' S[z] and Sq[z] 
are given by 



S[z] = S [z] + nx[z] 



dr z{t) z 



(19) 
(20) 



Let us call Qi{t) the ratio of path integrals appearing in 
Eq. JTHJ}, i.e. Q{t) = exp(-Ait 1/2 )Qi(i). The calculation 
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of the leading-order contribution to the asymptotic be- 
havior of Q\{t) can be obtained by restricting the sums- 
over-paths to those which start and end at z = 0, i.e. 
paths for which z[t) = z(Q) = 0. We can see this as 
follows. 

The expression for Qi(t) can be rewritten in the form 



i(*) = 



SI 



dxJ z z ^Z X Dz{t)eM-S[z]) 
I-oo dx^Zl Dz(t) exp(-SoH) 



(21) 



Let z c (x, t) and z co (x,t) be the trajectories that min- 
imize the functionals S[z] and Sq[z] respectively for 
boundary conditions z(0) = 0, z(t) = x, and write 
z(t) — z c (x, t) + z{t) in S[z] and z(r) = z co (x, t) + z(r) 
in S'ol-z]. Because both functionals are quadratic, we ob- 
tain immediately 



Qi(t) = 



I-oo dx exp{-S[z c (x, t)]} 
I-oc dx exp{-5 [z co (x,<)]} 

X Im=o D m^M-So[z})' 



(22) 



In the first factor in l|22|l . the actions S and So are 
both proportional to x 2 . For the 'free' action, So[z] — 
(1/4Da) L z 2 dr, the path z co is given by z co {x, r) = xr/t 
and the corresponding action is So^co] = x 2 /AD At- For 
the action S[z] one can show Q that S[z c ] oc x 2 /t^ with 
<f> > 1/4. The integrals over x in the numerator and 
denominator of the first factor in 1221) therefore yield just 
powers of t. By contrast, the second factor gives the 
exponential of a positive power of t. This second factor 
provides the leading correction to asymptopia, and will 
therefore be the focus of the subsequent development. To 
summarize, therefore, we will investigate the asymptotic 
behavior of the quantity 



i(*) 



fl^ZlD~z{t)eMS[S\) 



f^ZlD~z{t)eMS a [z\) 



(23) 



Since the trajectories z(t) vanish at both r = and r = t 
we can express them as a Fourier sine series: 

oo 

2(t) = a n sin(n7rr/t) (24) 



drz{r) sm(niTT/t). (25) 



We substitute the expansion (|2^|) into Eq. J22Jl, where 
the functionals S and S are given by l(2"U|) and pi[z] 
by l(T7|) . After the rescaling a n — ► 2[(D J 4i) 1 / 2 /im]a n we 
have: 



W)=ff[{^)e-^<-^ 

n—1 V 



(26) 



where the dimensionless coupling constant g is given by 

(27) 



D A Xit 1 ' 2 _ Ap D A 1/2 



7T 7 /2 £,1/2 



and the matrix elements A m „ are given by 



dy 



A mn = -1- / — 

n^Wo Vl-xJo y 1/2 {x~y) 3 / 2 

x[sin(m7ra;) — sin(mTry)] [sin(mrx) — sm(mry)]. 

(28) 

We see that the time dependence in the sub-leading 
contribution Qi(t) to Q(t) enters only through the 
coupling constant g. We note that the prefactor 
{Ap/Tr 7 l 2 ){D A /D)! 2 ) is small even for D A = D B = 1 and 
p = 1, where the 'pre-asymptotic regime' will be reached 
only for t 3> 7r 7 ~ 3 x 10 3 . At this timescale, the lead- 
ing asymptotic result is Q(t) ~ exp(— Xit 1 / 2 ) ~ 10~ 30 . 
This already gives strong hints as to why the asymptotic 
regime is so hard to reach. 

A full solution of the problem, at the quadratic level, 
requires the spectrum of the matrix A mn . At the moment 
we lack such an exact solution, so instead we devise a two- 
step approach to the problem: (i) we obtain an analytical 
lower bound for Q\(t); (ii) we calculate numerically the 
spectrum of the matrix A, truncated at a finite dimension 
N, and we make a finite-size scaling analysis for the sub- 
leading correction Q%(t,N) at finite N. 



C. A Lower Bound on Qi(t) 

Consider Eq. 126|l . By extracting the diagonal terms, 
(5/2) Anndnj from the quadratic form in the expo- 
nential, this equation can be rewritten as 



Qi(t) = n( 1+ ^«»r 1/2 (exp 



n=l 



(29) 

where the average is performed with the Gaussian weight 



exp 



IX^ 1 + 9Ann)a 2 t 



(30) 



Here we are using the diagonal part of the action, in the 
Fourier basis, as a trial action functional in a variational 
calculation. The variational property follows from the 
convexity inequality, (expx) > exp(x), which immedi- 
ately implies 



Qi(t) > Qu(t), 



where 



l (t) = exp 



i^ln(l+. gJ 4 n „) 



(31) 



(32) 
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is a lower bound for Qi(t). In the Appendix we show 
that, for large n, A nn ps ^/2n 2 n~ 3 / 2 . For large g, corre- 
sponding to large t, the sum on n is dominated by large 
values of n, of order n ~ <? 2 / 3 ~ i 1 / 3 , and the sum can be 
replaced by an integral to leading order: 

1 



hi(t) 



exp 
cxp(~Ag 2 / 3 ) 



-]Tln(l + .g\/27r 2 n- 3 / 2 ; 



where 

A = 2 1 / 3 tt 7 / 3 3- 1 / 2 . 
Using Ea.l|27|) for g, the final result takes the form 

l2 \ 1/3" 



Qu(t) w exp 



-4= ( 



(33) 
(34) 

(35) 



It is interesting to compare this result with the rigorous 
lower bound l)15[l. The latter is valid for all Da, while the 
new bound is valid, in principle, only for Da <C Db- In 
both expressions the leading correction to asymptopia is 
a term of order (p 2 Dt) 1 / 3 in the exponential, where D has 
the dimensions of a diffusion constant. In the bound H15(l 
D = D A , while in the new bound fJSJ L> = D 2 A /D B . For 
small enough Da/Db, therefore, the new bound is tighter 
than the old, while for large Da/D b the old bound is 
tighter. This confirms that the new bound cannot be 
valid in general, but only for small enough Da/D b . The 
two bounds cross when D A /D B = 3 9 / 2 tt 2 /128 ps 10.82, 
so the new bound is better for the case Da/D b = 1 used 
in simulations. Although the new bound is not strictly 
valid for Da = Db since the condition Da <C Db is not 
satisfied, the above considerations suggest that Da <C 
IODb may be sufficient in practice. 

Data for Q(t) were obtained using the algorithm of 
Mehra and Grassberger . The system parameters were 
D A = D B = D = 1/2 and p = 1/2. We plot, in Figure 
1, - h\Q(t)/(p 2 Dt) 1 / 2 against \n(p 2 Dt). The asymptotic 
value of this quantity is a\ = ps 2.257. The upper 

and lower bounds l (T4*|l and <f]~5f) are represented by the 
lower and upper dashed lines respectively (note the ap- 
parent interchange of bounds due to the minus sign in 
the definition of the ordinate). The new lower bound 
(|35|l is represented by the dotted line. It is clearly a sig- 
nificant improvement over the old bound, and appears to 
be converging towards the data at late times. 

In the following subsection we provide a finite- 
size scaling analysis which strongly suggests that the 
lower bound, Qu{t), gives the correct asymptotic form, 
Qi(t) ~ exp(— ci 1 / 3 ), although the constant c may differ 
from (i.e. be smaller than) that appearing in Eq. I|35|) . 



D. Exact Diagonalization for finite N 

One can show that the matrix elements A mn decay as a 
power law along the rows, so it is not obvious whether or 



— Numerical data 
— ■ Bounds 



Q 

a 


a 



L 



-2 







2 L 4 

ln(p"Dt) 

FIG. 1: Bounds on the quantity — InC 
tinuous line: numerical data from ref.|; 
D — Da = Db ~ 1/2; dashed lines: 
lower bounds; dotted line: new lower bound on Q(t) valid for 
D A < D B . 



)(t)/(p 2 Dt) 1/2 . Con- 
I, with p = 1/2 and 
rigorous upper and 



not the off-diagonal matrix elements make a significant 
contribution to the eigenvalues. These off-diagonal terms 
do not enter the calculation of the lower bound Qu, so it 
is clearly of interest to know whether they qualitatively 
change the result. 

As an analytical insight is missing so far, we study 
this problem numerically using the eigenvalues of finite 
matrices A^ n to compute Qi(g, N) of Eq. J2?J. The 
latter equation gives 



= [det(I + gA)} 



-1/2 



= exp | -- ^ ln(l + gX n ) j 



(36) 



where I is the unit matrix and A„ is the nth eigenvalue of 
A (where the eigenvalues are listed in decreasing order). 
Ideally, in a finite- N analysis we would simply truncate 
the sum at n = N. In practice, of course, we do not know 
the eigenvalues explicitly, so instead we truncate the ma- 
trix A to an N x N matrix and compute the eigenvalues 
of the truncated matrix. 

In the asymptotic regime we anticipate the finite-size 
scaling form 



HQi(9,N)] = N a f(g/N f3 ) . 



(37) 



The condition that Qi(g, N) has an N- independent large- 
TV limit implies that f(x) ~ x a ^ for x — * 0. 

It is instructive to first consider the diagonal ap- 
proximation to the matrix A, in which the off-diagonal 
terms are neglected. This provides the lower bound 
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Qu given by Eq. The corresponding finite- 

size quantity, Qu{g,N) is given by Qu(g,N) = 
exp[-(l/2) Y, n =i ln ( 1 + 9 Ann), leading to 

1 N 

- HQii(9, N)} = - ln (! + 9A nn ) ■ (38) 

1 n=l 

For large n, we show in the Appendix that A nn has the 
asymptotic form A nn ~ n~ 3 / 2 , so the finite-size scaling 
variable is g/N 3 / 2 , i.e. (3 — 3/2 in within the diag- 
onal approximation. In the limit g — > oo, N — > oo with 
g/N 3 / 2 held fixed, one can replace the sum over n by an 
integral and one readily recovers the scaling form 1|37|) 
with a = 1, (3 = 3/2. It is also easy to show that the 
scaling function f(x) in Eq. (|37|l has, within the diagonal 
approximation, the limiting forms fix) ~ x 2 ! 3 for 
as expected, while f(x) ~ ln(x) for i> 1. 

Our aim is to investigate whether the same general 
scaling structure (|37|l . with the same exponents a = 1 
and (3 = 3/2, holds when the full matrix A is used in- 
stead of just its diagonal part. If the matrix A is sim- 
ply truncated to an N x N matrix, however, and the 
eigenvalues determined, the corrections to finite-size scal- 
ing are found to be large. We attribute this to a large 
influence of the missing matrix elements A mn , with m 
or n greater than N, on the eigenvalues A„ when n 
gets close to N. Instead, therefore, we adopt a more 
refined approach in which the matrix is truncated to 
size 2N, but we use only the first N eigenvalues to 
compute Qi(g,N). Furthermore, instead of computing 

— ln<5i(.g, N) — (1/2) Xm=i m (l + 9^n), we compute the 
derivative 

H(g,N) = -2-fkiQxig, N) 
ag 

N 

= ^A„/(1+. 9 A„) , (39) 

n=l 

since we find that corrections to finite-iV scaling are 
smaller for this quantity. 

The finite-size scaling form for H(g, N) follows from 
its definition and Eq. I|37|l : 

H(g,N) = N a -^h(g/N^) , (40) 

where h(x) — 2df/dx. On the basis of the diago- 
nal approximation we expect a — 1, (3 = 3/2, and 
h{x) — ► (4A/3)x _1 / 3 for x <C 1, where A is given by 
Eq. 0, and, from Eq. h(x) -> x^ 1 for x » 1. We 
therefore plot N l / 2 H{g, N) against the scaling variable 

x = g/N 3 / 2 . 

Fig. 121 shows that a good scaling collapse is obtained 
with a = 1 and (3 — 3/2. This result implies that 

— hxQi oc g 2 ! 3 oc t 1 / 3 , i.e. the nondiagonal elements of 
the matrix A do not qualitatively change the scaling of 
the pre-asymptotic correction to the survival probabil- 
ity. The straight lines in the Figure are the asymptotes 
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FIG. 2: Scaling plot for H(g,N) = -2(d/dg)ln[Qi(g,N)] 
(N = 50,100,200,400,800,1600), with scaling variable x = 
g/N 3 ^ 2 . Continuous curves - full matrix; broken curves - 
diagonal matrix. The straight lines are the asymptotes for 
large and small scaling variable - see text. 



dicussed above, with slopes —1/3 and —1 for small and 
large scaling variable respectively. The main curve actu- 
ally contains both the full matrix (continuous lines) and 
the diagonal approximation (symbols), which lie almost 
on top of each other. 

A very small deviation from perfect scaling is evident 
in the data. In fact, a slightly better fit is obtained with 
(3 = 1.4 instead of 1.5. However, the same numerical scal- 
ing analysis applied to the diagonal matrix, where a = 1 
and (3 = 3/2 are exact, also yields, for the range of and 
g covered by Fig. (J2J, an apparent best fit with j3 « 1.4. 
We conclude that, for the full matrix, (3 = 3/2 is consis- 
tent with the data within the numerical precision. We 
conclude, finally, that the form of the leading correction 
to asymptopia given by the lower bound 133fl is correct. 
Even the amplitude seems to be given rather accurately 
by the lower bound. 



V. BEYOND THE QUADRATIC ACTION 

In the previous sections we have used the functional 
Hi[z] correct to first order in the ratio Da/Db of the 
diffusion constant of the particle to that of the traps. 
To this leading order, the parameter which controls the 
asymptotic behaviour is g oc Da/\TDb- We obtained 
the following asymptotic lower bound for the survival 
probability, Q(t), of the A-particle in d = 1: 

Q(t) > exphAxt 1 / 2 - Ag 2 ' 3 + ■■■} (41) 
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where g = (Da/Db)^ 1 ^ 2 /tt 3 and A is a constant. This 
suggests, as discussed in the preceding section, that at 
large Da/Db the new lower bound can cross the lower 
bound deduced from an "excluded volume" kind of argu- 
ment 0, y|, given by Eq. (JTSJ. 

To escape from this apparent contradiction, one has to 
consider higher-order (in Da/Db) terms in the action. 
To go to next order, we start from Eq. i|13[) . and expand 
the right-hand side up to order z 4 . After a tedious but 
straightforward calculation, we obtain a result for the 
term order z A in the form of a sum of (multiple) inte- 
grals. Inserting the Fourier expansion l|25l) and rescaling 
the expansion coefficients before to make them di- 

mensionless, we can cast the action in the dimensionless 
form 



S = S + gS! 



g^S 2 + 0([D A /D B } 2 
Db 



where 



So — 2 ^2 a » 

n 

S\ ~ ^ ^ A rnn Cl rn (l n 
m,n 

S2 ~ ^ A mn pqd rn CL n (lpC 



(42) 

(43) 
(44) 
(45) 



and the vertices A mn , A mnpq are pure numbers. Eq. I|42|) 
enables us to extend the calculation of Qi(t), Eq. I|26|) . to 
one higher order in Da/Db by replacing the exponential 
in Eq. $3% by exp(-S') with S given by Eq. g2J>- 

While we have not carried out this program explic- 
itly, the structure of the expansion l|42() shows that 
the quadratic approximation is valid for small Da/Db- 
What is the effect of the higher order terms? For the 
case Db = (where the traps are immobile) the asymp- 
totic behavior has, in d = 1, the known functional form 
exp(— ci 1 / 3 ), where c oc pD 1 ^ 3 . It is plausible to assume 
that higher-order corrections change the coefficient c but 
leave the exponent unaltered. Then, if the above series in 
Da/Db is convergent, we obtain the asymptotic behavior 
of Q(t), including the leading pre-asymptotic correction, 
in the form 



Q(t) ~ exp 



-^(p 2 D B t) 

\/7T 



1/2 



(46) 

The function F(x) has the limiting behavior F(x) — > 
Aix 1 / 3 for x -> 0, with A x < (32) 1 / 3 /\/3 from the 
lower bound i|35|) . while F(x) — > A2 for x — > 00, where 
A2 = 3/2 2 / 3 to match the known asymptotic behavior 
with immobile traps - the so-called Donsker-Varadhan 
problem [llj. This latter limit agrees with Eq. 115|) 
(where p$ = for immobile traps). The reason is sim- 
ple. In d = 1, only the two immobile traps immedi- 
ately to the left and right respectively of the A parti- 
cle play any role. Let the distances of these traps from 



the ^4-particle be x and y at t = 0. Then the parti- 
cle survives as long as it stays within the box of size 
(x + y) defined by these two traps, so its asymptotic 
survival probability, averaged over x and y, is given by 
Q(t) ^ P 2 Jo° C dx Jo°° dyexp[-p{x + y)-ir 2 D A t/(x + y) 2 }. 
Evaluating the integral by steepest descents for t — > 00 
one recovers Eq. (|15|l (with po — 0). 



VI. SUMMARY AND CONCLUSION 

In this paper we have addressed the issue of the lead- 
ing correction to asymptopia in the long-time dynamics 
of the trapping reaction, A + B — > B. The main re- 
sults are a new rigorous lower bound, valid in the limit 
Da <C Db, on the survival probability of the A-particle, 
and a numerical analysis, based on finite-size scaling, in- 
dicating that the new lower bound contains the correct 
analytic form for the correction to asymptopia, namely 
a i 1 / 3 correction to the leading i 1 / 2 term in the expo- 
nential. The new bound agrees well, at late times, with 
numerical data, as shown in Figure 1. We also presented 
a conjecture, Eq. I|46fl. for the form of the correction to 
asymptopia for general Da/Db, incorporating both ex- 
treme limits, Da <C Db and Da S> Db- 

Open questions to be addressed in future work include 
the generalization of the work presented here to higher 
spatial dimensions. Another interesting question, ad- 
dressed briefly in Q , is the scaling with time of the fluc- 
tuations of surviving trajectories. These can be studied 
starting from Eq. (|17|) , using the same Fourier decompo- 
sition of the trajectory as employed in this work. One 
finds that the fluctuations are subdiffusive as suggested 
in and observed in numerical simulations . Details 
of these calculations will be presented in a separate pub- 
lication. 



Acknowledgements 

We thank Richard Blythe for discussions. LA acknowl- 
edges support from the European Community Marie 
Curie Fellowship scheme under contract No. HPMF-CT- 
2002-01910. 



APPENDIX A: DIAGONAL ELEMENTS OF THE 
MATRIX A 

The matrix element A mn is given by the double integral 

/ -mnA t dX [ X dV 

x (sin(n7ra) — sin(n7r?/))(sin(m7ra;) — sin(m7r?/)) 

(Al) 

which by variable substitutions z — x—y, s = y and stan- 
dard manipulations of the trigonometric functions can be 
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reduced to a sum of one-dimensional integrals 

n(m — n) f 1 1 . nixz muz 
' dz ^tttt sin — — sin ■ 



2tt cos ■ 



xJo (|(m-n)(l-z)) 

ir(m + n) f 1 1 . nnz . muz 
' dz „ ,„ sin - sin - 



+ 27T COS ■ 



x Jo (^(m + n)(l - z)) 



(A2) 



where Jo (2) is the Bessel function of the first kind. The 
fact that I mn = Oifm + n = 2p + 1 can be obtained 
from the initial formula, Eq. (|A1|) , if one notices that the 
kernel is symmetric at reflection about the line y = 1 — x. 



For the diagonal elements, m — n, the leading large-n 
contribution comes from the first integral 



. ' 1 . 9 mr* 

277 1 dZ ^3j2 Sm — 



lo z 

/o 2 1/2 

v 2tt n ' 



1 



37i sm T 



(A3) 



correct to leading order for large n. It follows that 

A nn -> ^rr 2 n- 3 / 2 (A4) 

for n — > 00. 

The calculation of the asymptotic behaviour of A mn at 
large n with m fixed and large is more elaborate. Since it 
is not needed here, we will defer it to a future publication. 
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